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Abstract 

The authors present a relativistic and cross-section factorized framework for computing quasielas- 
tic A{p,pN) observables at intermediate and high energies. The model is based on the eikonal 



J2 ' approximation and can accomodate both optical potentials and the Glauber method for dealing 



with the initial- and final-state interactions (IFSI). At lower nucleon energies, the optical-potential 
philosophy is preferred, whereas at higher energies the Glauber method is more natural. This ver- 
satility in dealing with the IFSI allows one to describe A{p,pN) reactions in a wide energy range. 



. Most results presented here use optical potentials as this approach is argued to be the optimum 

m 

' choice for the kinematics of the experiments considered in the present paper. The properties of the 

C3 ' IFSI factor, a function wherein the entire effect of the IFSI is contained, are studied in detail. The 
predictions of the presented framework are compared with two kinematically different experiments. 

O I First, differential cross sections for quasielastic proton scattering at 1 GeV off ^-^C, ^^O, and '^'^Ca 

^ ' target nuclei are computed and compared to data from PNPI. Second, the formalism is applied to 



the analysis of a ^IIe(p, 2p) experiment at 250 MeV. The optical-potential calculations are found to 
be in good agreement with the data from both experiments, showing the reliability of the adopted 
model in a wide energy range. 
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FIG. 1: Schematic representation of the A{p,2p) reaction. The incoming proton undergoes "soft" 
initial-state interactions with the target before knocking out a bound proton through the primary 
high-momentum-transfer pp scattering. Both the scattered and the ejected proton suffer final- 
state interactions while leaving the nucleus. The scattered and the ejected proton are detected in 
coincidence, while the residual nucleus remains unobserved. 



I. INTRODUCTION 



Quasielastic nucleoli knockout reactions have been extensively investigated with the aim 
of obtaining precise information on nuclear structure. The present work focuses on exclusive 
proton-induced A{p,2p) and A{p,pn) processes, whereby the residual A-1 nucleus is left in 
the discrete part of its energy spectrum. A sketch of the A{p, 2p) reaction is given in Fig. ^ 
For a historical overview of the research into proton-induced nucleon emission off nuclear 
targets the ..eade.,s..efe„.ed to Refs.QBa. 

In quasielastic A{p,pN) scattering the projectile is elastically scattered from a single 
bound nucleon in the target nucleus, resulting in the struck nucleon being knocked out of 
the target nucleus. This relatively simple reaction mechanism of one "hard" nucleon-nucleon 
collision is obscured by the "soft" initial- and final-state interactions (IFSI) of the incident 
and two outgoing nucleons with the nuclear medium. Consequently, every model for A{p, 2p) 
and A{p,pn) reactions has two issues to address : first, the description of the "hard" wide- 
angle scattering that leads to the ejection of the struck nucleon and second, the distorting 
mechanisms of the "soft" small-angle IFSI. 

Concerning the treatment of the "hard" A^A^ scattering part, essentially two methods 
exist. A so-called "cross-section factorized" approximation adopted so that the 

nucleon-nucleon scattering cross section enters as a multiplicative factor in the differential 
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A{p, pN) cross section. Some results of exclusive Alp^p) measurements interpreted with this 



BQBBQ. 



cross-section factorized form can be found in Refs. IM l6l LZl |8|, |9|, lldl • The inclusion of spin- 
dependence in the description of the IFSI, however, breaks this factorization scheme. In that 
situation an alternative technique can be used : the amplitude factorized form of the cross 
section Q]. In this approach, the two-body A^A^ interaction can be approximated by the 
interpolation of phase shifts from free elastic A^A^ scattering. Various phenomenological 
forms to fit the amplitudes have also been used in the past. Traditionally, the nucleon- 



nucleon scattering matrix has been parametrized in terms of five Lorentz invariants 



a method usually dubbed as the lAl model or the SPVAT (scalar, pseudoscalar, vector, 
axial vector, tensor) form of the A^A^ scattering matrix. Differential cross section calculations 
adopting these five-term representations have been reported in Refs 

Q H Q, Q, Q- 

It should be noted, however, that although the SPVAT form gives reasonable predictions 
of A{p,pN) observables, it is, in principle, not correct, as a five-term parametrization of 
the relativistic A^A^ scattering matrix is inherently ambiguous [2^. Tjon and Wallace [2^ 
have shown that a complete expansion of the A^A^ scattering matrix (commonly called the 
IA2 model) contains 44 independent invariant amplitudes. To date, the only calculations 
employing this general Lorentz invariant representation have been performed in the context 
of the relativistic plane wave impulse approximation (RPWIA), i.e., a model which ignores 
all IFSI mechanisms 123^ . 

The IFSI effects are typically computed by means of the distorted wave impulse ap- 



proximation (DWIA) theoretical framework 




Generally, in a DWIA approach 



the scattering wave functions of the incoming and two outgoing nucleons are generated by 
solving the Schrodinger or Dirac equation with complex optical potentials. Parametriza- 
tions for these optical potentials are usually not gained from basic grounds, but are ob- 



tained by fittin g e 



eter sets 



astic nucleon-nucleus scattering data. Several optical-potential param- 



27L I28L I29L l3flL I31l | have been used in the description of quasifree pro- 
ton scattering, off nuclei. In the past, both nonrelativistic and relativistic DWIA versions 
0, y, fl J ID, Q 12, H Q, H H 3, Q, 3 have proven successful in predicting 
A{p,pN) cross sections over a wide energy range (76-600 MeV) and for a whole scope of 
target nuclei. 

Most of the available calculations for the exclusive A{p,pN) process addressed incident 
proton kinetic energies of a few hundred MeV. In this work we aim at extending the for- 
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malism to scattering in the GeV energy regime. As a matter of fact, the majority of DWIA 
frameworks rely on partial-wave expansions of the exact solution to the scattering problem, 
an approach which becomes increasirigly cumbersome at higher energies. In this energy 
range the eikonal approximation js^, IssI' which belongs to the class of high-energy semi- 
classical methods, offers a valid alternative for describing the IFSI. Nonrelativistic eikonal 
studies of the A{p, 2p) reaction used in combination with optical potentials can be found in 
Refs. IllflQ- 

In this paper we propose a relativistic and cross-section factorized formalism based on the 
eikonal approximation for computing exclusive A{p,pN) cross sections at incident proton 
energies in the few hundred MeV to GeV range. The eikonal formalism is implemented 
relativistically in combination with optical potentials js^, as well as with Glauber theory 



3S 



40 



4l| , which is a multiple-scattering extension of the eikonal approximation. The two 



frameworks only differ in the way they treat the IFSI. The main focus will be on the optical- 
potential approach, as this method turns out to be the most suitable for the description of 
the IFSI for the kinematical settings discussed in this work. 

The paper is organized as follows. In Sees. Ill Al and IIIBI the factorized cross section 
is derived in the RPWIA formalism. Thereafter, the different methods to deal with the 
IFSI are developed in Sec. Ill Gl Sec. IIIII is devoted to a presentation of optical-potential 
and Glauber results for the IFSI factor. This is a function which accounts for all IFSI 
effects when computing the A{p,pN) observables. The optical-potential predictions of our 
model are compared with cross section data that have been collected at PNPI and TRIUMF 
in Sec. II VI First, we present our calculations for the ^^G, ^^O, and ^'^Ga(p, 2p) and {p,pn) 
PNPI data for 1 GeV incoming proton energies 42|. Second, the cross sections for ^He(p, 2p) 
scattering at 250 MeV are compared to the TRIUMF data by van Oers et al. j^. Finally, 
Sec. IVl states our conclusions. 



II. A{p,2p) FORMALISM 



In this section the formalism for the description of A{p, 2p) reactions is outlined. The 
generalization to AJmpn) reactions is straightforward. We conform to the conventions of 
Bjorken and Drell 43[ for the 7 matrices and the Dirac spinors, and take h = c = 1. 
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A. A{p, 2p) differential cross section and matrix element 



The four-momenta of the incident and scattered proton are denoted as Pi{Epi,pi) 
and Ki{Eki, ki). The proton momenta pi and ki define the scattering plane. The four- 
momentum transfer is given by {u, q) = = Pi^—K'^ = Kj^_-^+K2—K^, where K^{Ea, kji), 
K'^_^{Ea-i, kA-i), and (-E'fc2, ^2) are the four-momenta of the target nucleus, residual 



nucleus, and the ejected proton. The standard convention Q" = —q^q^ 
followed for the four-momentum transfer. 

In the laboratory frame, the fivefold differential cross section can be written as 



cu^ > is 



a 



dEkidVtidVt'2 



kik2 .-1 v^l 



(P,2p) 
fi 



if 



(1) 



Here, Al^^'^^^ is the invariant matrix element which reflects the transition between the initial 
and final states. The hadronic recoil factor is given by 



fr 



Ea- 



Ea 



1 



E^ 



k2 



Ea^ 



q- k2 

^2 



1 



Ujk2 - qEk2 cos Ok^q 



MAk2 



(2) 



with the energy transfer u = Epi — E^i = Ea-i + Ek2 — Ea, the three-momentum transfer 
q = Pi — ki = kA-i + k2 — kA, and Ok^q the angle between k2 and q. 
The A{p, 2p) matrix element is given by 



(i^fm,!^, K2^m,2/, A-1 {K>X_^, Mr) | O^^) |p/^ m,i,, A (J^^, 0+, g.s.)) 



(3) 



(4) 



where 

j<i=o 

is the unknown two-body operator describing the high-momentum transfer "hard" pp 
scattering, \A (i^^, 0^, g.s.)) the ground state of the even-even target nucleus and 
|A-1 JrMr)'^ the discrete state in which the residual nucleus is left. In coordinate 

space the matrix element takes on the form 



M 



fi 



dfo I dfi I df2... I rff^4 ( ^^''T"'-^'''''''""^-^ (ro,n,r2. 



fA))'d^'^ 



xn+?"''' (ro,n,f2,...,f^) 



(5) 



For the sake of brevity of the notations, only the spatial coordinates are explicitly written. 



B. Relativistic plane wave impulse approximation 



In this section, the A{p, 2p) matrix element of Eq. will be analyzed in the RPWIA. 
In this approach, only one hard collision between the projectile and a bound nucleon is 
assumed to occur, knocking the bound nucleon out of the target nucleus. The modelling 
of the "soft" IFSI processes, which affect both the incoming and outgoing protons, will be 
considered in Sec. Ill Ul 

In evaluating the A{p, 2p) matrix element of Eq. (0), a mean-field approximation for the 
nuclear wave functions is adopted. We also assume factorization between the "hard" NN 
coupling and the nuclear dynamics. For reasons of conciseness, the forthcoming derivations 
are explained for the A = 3 case. The generalization to arbitrary mass number A is rather 
straightforward. 

The antisymmetrized (A+l)-body wave function in the initial state is of the Slater de- 
terminant form 



ViA + iy. 



(ro) (par (ro) (fo) (pas (ro) 

0Pim,n in) (par (fi) (pa^ {n) (pa^ {ti) 

in) (Pai in) (Pa2 (^2) (Pa-i i^) 

(^3) <Pai (^3) <Pa2 (^3) <Pa3 (^3) 



(6) 



Details on the bound-state single-particle wave functions 0q,- (r*, a) entering this mean-field 
(A+l)-body wave function can be found in Appendix^ The wave function of the incoming 
proton is given by a relativistic plane wave 



kms 



[r = 



E + M 
2M 



Ak-f 



u{k, rris 



(7) 



The (A+l)-body wave function in the final state reads 



A+1 



'ro,fi,f2,f3) 



knusif 


(ro) 


^k2ms2f 


klTUslf 


(ri) 


^k2ms2f 




in) 


^k2ms2f 


kinisif 


in) 


^k2ms2f 



in) (Pa2 in) (pas in) 



(8) 



Relative to the target nucleus ground state written in Eq. (jH)), the wave function of Eq. (jH} 
refers to the situation whereby the struck proton resides in a state "ai" , leaving the residual 
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A-1 nucleus as a hole state in that particular single-particle level. The outgoing protons are 
represented by relativistic plane waves. 

Since both the initial and the final wave functions are fully antisymmetrized, one can 
choose the operator O^"^^ to act on two particular coordinates (ro and fi). Without any loss 
of generality, the A{p, 2p) matrix element of Eq. can be written as 

^ E E E E 

X O (fo, fi) (po (fo) (pp (ri) (pq (ra) (pr (fs) , (9) 
with eijki the Levi-Civita tensor. In the RPWIA, 

= 5mg j dro j dfi j (if20|.(fo)0[(ri)O(fo,ri)0o(fo)0p(fi) |0g(f2)|^ . (10) 
Inserting this expression in Eq. Q one obtains 

X E EE 

X (^klmn^opmn 

0l(^"*o)0j'(ri)|0„^(rl)H0„(rl)|'O(rl,,ri)</.o(ro)</.p(ri) . (11) 

There are {A — 1)\ possible choices (permutations) for the indices m, n, . . . , all giving the 
same contribution to the matrix element. Accordingly, the above expression can be rewritten 

as 

Mff'^ = I J dfo J dn J df, J df, J2 E 

k,l€^kim.sif ,k2ms2f} o,pe{pimsH,«i} 

4>i (^o) 4>\ [ri] (fa)!^ 1003 (fg)!^ O (ro, fi) (po (fo) (pp (fi) (12) 

Because the bound-state wave functions are normalized to unity (J df\<pa{r)\^ = 1) and 
O (fo, fi) = O (fi, fo), the matrix element can be further simplified to 

>i(f = / dro [ dn Ul (fo) 4 (^"i) - 4 (^"o) 4 (^"i)) 

/* J J \ kirusif ^ ' k2ms2f ^ ' k2ms2f ^ ' ^ kmisif ^ ' J 

X O (fo, fi) 0pl„i,i, (fo) 0ai (fl) , (13) 



including a direct and an exchange term. 

Substitution of the general form of the scattering operator O {fo,fi) = J j0^e'^^'^'^^~^°''F, 
with F the NN scattering amplitude in momentum space, in the above expression (fT!?|l leads 

to 

= v) {ki,msif)u\k2,ms2j) F u{pi,msii)(j)ai (Pm) (14) 

for the direct term and a similar expression for the exchange term. Here (pa (p) is the rela- 
tivistic wave function for the bound nucleon in momentum space (for details see Appendix IX| 
and Pm = ki + k2 — pi is the missing momentum. In order to arrive at a cross-section factor- 
ized expression for Eq. the quasielastic off-shell proton-proton scattering matrix element 
will be related to the free on-shell proton-proton cross section. For this purpose, we insert 
the completeness relation 

[u{Pm, S)u{prn, s) - v{pm, s)v{pm, s)] = 1 (15) 

s 

in 

Mff^^ = U^{ki,mslf)u^{k2,ms2f) F u{pi,msli)(l)aAPrn) - (kiTJlslf ^ k2ms2f^ , (16) 

and obtain the following expression for the matrix element 

s 

— negative-energy projection term . (17) 
Here, A^j^ is the matrix element for free pp scattering 

- (jtinisif ^ k2ms2f^ ■ (18) 

Factorization breaks down, even when IFSI are disregarded, owing to the negative-energy 
projection term. To recover factorization, the negative-energy projection term is neglected 
in the remainder of this work 

J<" « E (-^^S™...,,™,,,,™., s(p„, (p„.) . (19) 

S 
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Using the expression of the relativistic bound-nucleon wave function in momentum space 
given in Appendix |X1 the u (pa contraction in Eq. (fT^ reduces to 



2M„ 



(20) 



where E = ^"p^ + M^, xig3^Km indicates the spin projection of the spin spherical harmonic 
3^Km(^p) on a spin state Xi^! ^^^d the radial function in momentum space a^K is given by 



V 



E + M, 



(21) 



with and /„k the Bessel transforms of the standard upper and lower radial functions 
of the bound-nucleon wave function in coordinate space (see Appendix [X] for details) and 
Sk — f^l I'^l • 

Upon squaring Eq. (fT^ . the pp and nuclear bound-state parts get coupled by the sum- 
mation over the intermediate spins s and s' : 

s,s' 

X {u{pm,s)(j)aAPm)y u{Pm,s')(l)aAPm) ■ (22) 



M 



(P,2p) 
fi 



After summation over m, the struck nucleon's generalized angular momentum quantum 
number, the square of u{pm, s)(f)ai {Pm) yields a Sgs', i.e., becomes diagonal in s. Thereby, 
use is made of the following identity 

2j + l 



in 



(23) 



This leads to the decoupling between the pp scattering and the bound-state part in the 
matrix element : 



I'-'^nft (Pm) I 



X 



with 



ftriK (p) 



lE + Mp 
2M„ 



(ynnip) ■ 



(24) 



(25) 
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The last factor in Eq. ()24p can be related to the free pp scattering center-of-mass cross 
section 

("^^ =^lV Y\{Mn ' (26) 

SO that the RPWIA differential A{p, 2p) cross section of Eq. (Q) can be written in the cross- 
section factorized form 

RPWIA ^MA-ihk2 2j + l ,, M2(da^\ ..^^ 
\dE,,dQ,dQ2) ^M.MaPi^''^' 47r A.m. ' 

Here, s is the Mandelstam variable for the pp scattering, not to be confused with the inter- 
mediate spin from the preceding Eqs. In the numerical calculations which will 
be presented in the forthcoming sections, the free proton-proton cross section {^^)^^ is 



obtained from the SAID code 



C. Treatment of the IFSI 

It is well known that the factorized RPWIA result of Eq. ^I7\i adopts an oversimplified 
description of the reaction mechanism. The momentum distribution {Pm)f, which 

represents the probability of finding a proton in the target nucleus with missing momentum 
Pm, will be modified by the scatterings of the incoming and outgoing protons in the nucleus. 
Therefore it is necessary to incorporate the effects of these IFSI in the model. 

First, in Sec. Ill C 11 the differential A{p, 2p) cross section is written in a factorized form 
taking IFSI effects into account. Next, the relativistic eikonal methods used for dealing 
with the IFSI effects in this work, are discussed in depth. Two methods will be used. 
The relativistic optical model eikonal approximation (ROMEA) is the subject of Sec. Ill C 2| 
whereas the relativistic multiple-scattering Glauber approximation (RMSGA) is discussed 

in Sec. mm 



1. Factorization assumption and the distorted momentum distribution 

In both versions of the relativistic eikonal framework for A{p,pN) reactions presented 
here (ROMEA and RMSGA), the antisymmetrized initial- and final-state (A+l)-body wave 
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functions, 



'Spi (ro, f2, . . . , ta) e'^^""" u{pi, msu) (n, rg, . . . , f^) 



(28) 



and 



A+1 



ro,ri, . . . ,rAj 



A 



9"^ 



.ta e 



iki -fo 



u{ki,msif) 



Ta) e''=^-^^M(^2,m,2f) 



(29) 



differ from their respective RPWIA expressions of Eqs. © and (jH)) through the presence of 
the operators Spi, Ski, and Sk2- These define the accumulated effect of all interactions that 
the incoming and emerging protons undergo in their way into and out of the target nucleus. 

Since the IFSI violate factorization, some additional approximations are in order. First, 
only central IFSI are considered, i.e., spin-orbit contributions are omitted. Further, the 
zero-range approximation is adopted for the "hard" A^A^ interaction, allowing one to replace 
the coordinates of the two interacting protons (ro and fi) by one single collision point in the 
distorting functions Spi, Ski, and 5^2 • This leads to the distorted momentum-space wave 
function 

K (P-) = / dfe-'P^-Ua, (r) SjFsi (r) , (30) 
similar to Eq. ()A4|) . but with the additional IFSI factor 



S 



iFSiyr) 



dfo ... / df 



-A \Va2 (^2)1^ • • • (^vi)l^ Ski {r, ra, . . . , f^) 



X Sk2 (r, f2, . . . , Ta) Spi (f, f2,...,fA 



(31) 



accounting for the soft IFSI effects. 

Now, along the lines of 4^, it is natural to define a distorted wave amplitude 



D 



(32) 



so that the distorted momentum distribution is given by the square of this amplitude, 



(33) 



This distorted momentum distribution has the following properties. First, it takes into 
account the distortions for the incoming and outgoing protons. Second, it reduces to the 
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plane wave momentum distribution |a„K {Pm)f in the plane wave limit when assuming 
that (pai {Pm) satisfies the relation 



E + Mp 



(34) 



between the upper and lower components. 

Using the ansatz (j33|) for the distorted momentum distribution, the differential A{p, 2p) 
cross section can be cast in the form 

It differs from the RPWIA expression through the introduction of a "distorted" mo- 
mentum distribution p^. 



2. Relativistic optical model eikonal approximation 



As shown for example in Refs. 



46 



3, 



in the relativistic eikonal limit the scattering wave 



function of a nucleon with energy E = ^/k^ + and spin state subject to a scalar 

(Vg) and a vector potential {V^) takes on the form 



E + M 
2M 



1 



E+M+Vs-V^ 



a ■ p 



Jk-r iS{r) 

^ ^ Aims ' 



(36) 



where the eikonal phase S{h, z) reads 



M 

iS{h, z) = -I— 



dz' 



V,{h, z ') + Ko(6, z')[a-{h^K)- iKz '] 



(37) 



with r = (6, z) and the average momentum K pointing along the 2;-axis. The central and 
spin-orbit potentials Vc and Vso in the above expression are determined by Vg and K and their 
derivatives. In general, a fraction of the strength from the incident beam is removed from 
the elastic channel into the inelastic ones. These inelasticities are commonly implemented 
by means of the imaginary part of the optical potential. 

In evaluating the IFSI effects, three approximations are introduced. First, the dynamical 
enhancement of the lower component of the scattering wave function (jHE)), which is due 
to the combination of the scalar and vector potentials, is neglected. Second, the impulse 
operator p is replaced by the asymptotic momentum k of the nucleon. As mentioned before, 
the spin-orbit potential is also omitted. 
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As a result, the effects of the interactions of the incoming and outgoing protons with 
the residual nucleus are implemented in the distorted momentum-space wave function of 
Eq. (pUj) through the following phase factors 

Sp, (f) = e-'^ ""'^'"^ , (38a) 



ki [r 



(38b) 



5fc2(r) = e '^^^'2 \ (38c) 

with the 2;- axes of the different coordinate systems lying along the trajectories of the re- 
spective particles {z along the direction of the incoming proton pi, z' along the trajectory 
of the scattered proton /ci, and z" along the path of the ejected nucleon ^2) and (6p^,2;pJ, 
{bki,Zki), and {bk2,Zk2) the coordinates of the collision point r in the respective coordinate 
systems. The geometry of the scattering process is illustrated in Fig. |21 The integration 
limits guarantee that the incoming proton only undergoes ISI up to the point where the 
"hard" A^A^ collision occurs and the outgoing protons are only subject to FSI after this 
"hard" collision. 

It is worth remarking that the eikonal IFSI operators of Eq. ()38|) are one-body operators, 
i.e., they do not depend on the coordinates {r2,'r3, . . . ,rA) of the residual nucleons. The 
normalization of the bound-state wave functions simplifies the IFSI factor (jHTj) considerably 
to (5/^5/ (r) = Ski (r) Sk2 (r) Spi (r) in the ROMEA case. 

In the numerical calculations, we have employed the global S — V parametrizations of 
Cooper et al. and the optical potential of van Oers et al. to describe the PNPI and 
TRIUMF data, respectively. Hereafter, the A{p, 2p) calculations which adopt Eq. (jHHj) as a 
starting basis are labeled the relativistic optical model eikonal approximation (ROMEA). 



3. Relativistic multiple- scattering Glauber approximation 

In the ROMEA approach, all the IFSI effects are parametrized in terms of mean-field 
like optical potentials, i.e., the IFSI are seen as a scattering of the nucleon with the residual 
nucleus as a whole. As the energy increases, shorter distances are probed and the scattering 
with the individual nucleons becomes more relevant. For proton kinetic energies Tp > 1 GeV, 
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FIG. 2: Geometry of the scattering process. The vectors bp-^ , b^i , and bf^^ the impact parameters 
for each of the three paths for a coUision occurring at f. Zp-^ , , and are the z coordinates of 
the colhsion point in the respective coordinate systems. 9i and 62 are the angles of the outgoing 
nucleons relative to the incoming proton direction. 



the highly inelastic and diffractive character of the und erlying elementary proton-nucleon 
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4ll | more natural. This method 



scattering cross sections makes the Glauber approach [39 , 
reestablishes the link between proton-nucleus interactions and the elementary proton-proton 
and proton-neutron scattering. It essentially relies on the eikonal, or, equivalently, the 
small-angle approximation and the assumption of consecutive cumulative scatterings of a 
fast nucleon on a composite target containing "frozen" point scatterers (nucleons). 

In relativistic Glauber theory, the scattering wave function of a nucleon with energy 



E = ^/k^ + and spin state II^ti^) reads 

■,(+) _ 
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k,ms 



E + M 



2M 



S 



1 



1 



E+M 



a ■ p 



(39) 



The operator S implements the subsequent elastic or "mildly inelastic" collisions of the fast 
nucleon with the "frozen" spectator nucleons 

A 



S (r,f2,f3, . . . ,r^) 



ll[i-T(b-b,)e{z-z,) 

i=2 



(40) 



where 9 (z 



ensures that the nucleon only interacts with other nucleons if they are 



localized in its forward propagation path. Given the diffractive nature of pN collisions at 
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GeV energies, the profile function T ( 6 ) for central elastic pN scattering is parametrized in 
a functional form of the type 



At lower energies that part of the profile function proportional to ep^v is non-Gaussian and 



makes significant contributions to nuc 



ear scattering. Rather than Eq. (PT|) a parametrization 



in terms of the Arndt A^A^ phases ^| is appropriate at lower energies. For the calculations 
presented here, which address higher energies, the Gaussian-like real part of F {h^ is the 
dominant contributor and the use of Eq. (PT|) is justified. The parameters in Eq. (PT|) can 
be determined directly from elementary nucleon-nucleon scattering experiments and include 
the total pN cross sections 0"*'^^, the slope parameters j3pN, and the ratios of the real to 
the .„.agina.y part of the scattering a,npUtude .... We obtained these Glaube„ete. 
through interpolation of the data base available from the Particle Data Group j50|] (for more 
details, see Ref. As in the ROMEA framework, only the central spin-independent 

contribution is retained and the impulse operator is replaced by the nucleon momentum. 
The Glauber operators in Eq. take the following form 

A 

5pi(^^2,r3,...,rA) = n[l-r(^-^)^(^-%)] > (42a) 

i=2 

A 

Ski (f, f2, fa, . . . , f^) = J] [l - F (6 ' - ') e {z^ '-z')\ , (42b) 

A 

Sk2 (f, rl, fs, . . . , ta) = n - r - h ") ^ " - ^ ")] ' (42c) 

where f denotes the collision point and (7^2, fs, • • • , ^a) are the positions of the frozen specta- 
tor protons and neutrons in the target. The (6, 2;), (6', 2;'), and (6", z") coordinate systems 
are defined as in the previous section. The step functions make sure that the incoming 
proton can only interact with the spectator nucleons which it encounters before the "hard" 
collision and the outgoing protons can only interact with the spectator nucleons which they 
find in their forward propagation paths. 

Contrary to ROMEA, the Glauber IFSI operators of Eq. ()42|) are genuine A-body oper- 
ators, so the integration over the coordinates of the spectator nucleons in Eq. (|HT|) has to 
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be carried out explicitly. This makes the numerical evaluation of the Glauber IFSI factor 
very challenging. Standard numerical integration techniques were adopted to evaluate the 
IFSI factor and no additional approximations, such as the commonly used thickness-function 
approximation, were introduced. 

Henceforth, we refer to calculations based on Eq. (|42p as the relativistic multiple- 
scattering Glauber approximation (RMSGA). 

III. NUMERICAL RESULTS FOR THE IFSI FACTOR 

In this section, results for the IFSI factor (jHT|) are given for the knockout of nucleons 
from the Fermi level in ^^C, ^^O, and ^°Ca, at an incident energy Tpi = 1 GeV and a 
scattered proton kinetic energy T^i = 870 MeV. Thereby, we adopt coplanar scattering angles 
{01,62) = (13.4°, 67°) on opposite sides of the incident beam, i.e., kinematics corresponding 
with the PNPI experiment of Ref. [4^. All IFSI effects are included in the IFSI factor 
SiFsii"^- Note that in the absence of initial- and final-state interactions the real part of the 
IFSI factor equals one, whereas the imaginary part vanishes identically. 

The A{p, 2p) IFSI factor is a function of three independent variables (r, 6, (p). The z-axis 
is chosen along the direction of the incoming beam pi, the y-axis lies along pi x ki and the 
X-axis lies in the scattering plane defined by the proton momenta pi and ki. 9 and denote 
the polar and azimuthal angles with respect to the z-axis and the x-axis, respectively. The 
radial coordinate r represents the distance relative to the center of the target nucleus. 

A. 9 dependence 

To gain a better insight into the dependence of the IFSI factor on r, 9, and 0, we 
calculated the contribution of the three distorting functions iSpi, iS^i, and 8^2 to the IFSI 
factor. In Figs. El and |3] results are displayed for the computed real and imaginary part of 
<SiFsi{i^, 0,4' = 0), for proton emission from the Fermi level in ^^C. The results were computed 
within the ROMEA framework, using the EDADl optical-potential parametrization of j^^. 

The 9 dependence can be interpreted as follows. For a given r, the distance that the 
incoming proton travels through the target nucleus before colliding "hard" with a target 
nucleon, decreases with increasing angle 9. As a consequence, small values of 9 induce 
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FIG. 3: (Color online) The radial and polar-angle dependence of the real part of the IFSI factor 
<SiFSi in the scattering plane {(f) = 0°) for proton knockout from the Fermi level in ^^C. The upper 
left panel is the contribution from the impinging proton (Spi), while the upper right panel shows 
the effect of the FSI of the scattered proton {Ski)- In the bottom left figure, the effect of the FSI 
of the ejected proton (5^2) is presented and the bottom right figure shows the complete IFSI factor 
{Ski Sk2 Spi). The kinematics was Tpi = 1 GeV, Tki = 870 MeV, Oi = 13.4°, and 6*2 = 67°. 
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FIG. 4: (Color online) As in Fig. but now for the imaginary part of the IFSI factor. 



the largest ISI. For the FSI of the scattered proton, the opposite holds true and 6 = 180° 
for a large r value corresponds to an event whereby the "hard" collision transpires at the 
outskirts of the nucleus and the scattered proton has to travel through the whole nucleus 
before it becomes asymptotically free, thus giving rise to the smallest (largest) values for the 
real (imaginary) part of the IFSI factor. Unlike the scattered proton, which moves almost 
coUinear to the 2;-axis, the ejected nucleon leaves the nucleus under a large scattering angle 
6*2- Hence, the FSI are minimal for 6 close to 0° or 180° and maximal for 6 around 180° — 62. 
Finally, the 9 dependence of the complete IFSI factor is the result of the interplay between 
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FIG. 5: The radial dependence of the real part of the IFSI factor Sipsi in the scattering plane 
{(j) = 0°) for neutron knockout from the Fermi level in ^'^Ca. The upper left, upper right, and 
bottom left panels display Spi, Ski, and 5^2, respectively, while the complete IFSI factor is shown 
in the bottom right picture. 



the three distorting effects, with the strongest scattering and absorption observed at 6 close 
to 0°, 180° - 02, and 180°. 

B. r dependence 

Figure m displays the real part of the IFSI factor as a function of r at various values of 6. 
The ROMEA calculations were performed for the same kinematics as in Figs. El and EJ and 
use the EDAI optical-potential fit of 0]. 

Turning first to the upper left panel, it can be inferred from the picture that the ISI 
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effects increase with growing r for ^ = 0°. For 6 = 45° and increasing r, initially, the 
growing distance the proton has to travel through the nucleus leads to a decrease of the real 
part of Spi. This is followed by an increase for larger r up to Spi = 1. This reduction in ISI 
effects with increasing r is brought about by the incoming proton's path through the nucleus 
moving away from the nuclear interior and closer to the less dense nuclear surface. The other 
curves of the upper left figure illustrate a general trend for 90° < 9 < 180° : as r increases, 
the real part of the IFSI factor grows correspondingly. As can be appreciated from Fig. El 
as well as from the previous figures, the global behavior of the Ski factor describing the 
scattered proton's FSI can be related to that of the ISI factor Spi through the substitution 
9 — s> 180° — 9. This approximate symmetry can be attributed to the small scattering angle 
9i, i.e., the scattered proton leaves the nucleus almost parallel to the incoming proton's 
direction. In the bottom left panel, the additional curve {9 = 115°, i.e., close to 180° — ^2) 
represents the situation of maximal FSI of the ejected nucleon. For this 9 value, the path 
of the ejected nucleon passes through the center of the nucleus and the distance travelled 
through the nucleus increases with r. Accordingly, the real part of Sk2 is a monotonously 
decreasing function of r. The other extreme is the ^ = 0° case where increasing r means 
less FSI. For the other 9 values, the absorption reaches its maximum for some intermediate 
r value. Again, the combination of Spi, Ski, and Sk2 determines the total IFSI factor with 
the strongest attenuation predicted in the nuclear interior. 

C. 4> dependence 

The dependence of the IFSI factor on the azimuthal angle of the collision point is quite 
straightforward. One representative result is displayed in Fig. IHl Here, cos0 > (cos0 < 0) 
refers to a situation where the "hard" A^A^ collision occurs in the upper (lower) hemisphere 
with respect to the yz-pla.ne (see Fig. El for a collision point located in the upper hemi- 
sphere). Due to the cylindrical symmetry about the z-axis, the factor describing the ISI of 
the incoming proton is independent of 0. Regarding the scattered proton, we observe the 
least FSI in the upper hemisphere, since the proton then avoids passing through the highly 
absorbing nuclear interior. For the ejected nucleon, the contrary applies and the strongest 
FSI effects are found for = 0°. As the xz-plane is defined as the scattering plane, the IFSI 
factor possesses the symmetry Sipsii^, 9, 2ti — (p) = Sipsii^, 9, (p) for coplanar scattering. 
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FIG. 6: (Color online) The polar- and azimuthal-angle dependence of the real part of SjFsiif = 
3 fni,9,(j)) for proton knockout from the Fermi level in ^^O. Kinematics as in Fig. |31 ROMEA 
calculation with the EDAD2 optical potential [s]]. As in the previous figures, the upper left, upper 
right, and bottom left panels represent the effect of the ISI of the incoming proton, the FSI of the 
scattered proton, and the FSI of the struck nucleon, respectively, whereas the bottom right figure 
displays the complete IFSI factor. 
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D. Level and A dependence 



Results for the emission from levels other than the Fermi level have not been plotted here, 
as it turns out that the IFSI factors hardly depend on the single-particle level in which the 
struck nucleon resides. The peculiar spatial characteristics of the different single-particle 
orbits have an impact on the observables, though. Indeed, the distorted momentum-space 
wave function (p^^ of Eq. (jHUI) is determined by the values of the IFSI factor folded with a 
relativistic bound-state wave function (p^i (^)- As the particles experience less IFSI close to 
the nuclear surface, one obtains a stronger reduction of the quasifree cross section for nucleon 
knockout from a level which has a larger fraction of its density in the nuclear interior. This 
will become apparent in Fig. [7| but even more so in Sec. IIV Al 

Figure [7| shows a function 6 (r) which represents the contribution of the nuclear region 
with radial coordinate r to the differential cross section. The procedure to calculate this 

n 

function is similar to the method exposed in Ref. 51| and is developed in Appendix^ Com- 
parison of the upper and lower panels illustrates that IFSI mechanisms make the A{p, 2p) 
cross sections reflect surface mechanisms, unlike the A{e, e'p) reaction where the weakly in- 
teracting electron probes the entire nuclear volume and only the outgoing proton interacts 
with the residual system. Apart from the shift to higher r, the IFSI brings about a strong 
reduction in the magnitude of the cross sections, whereby the Fermi level is least affected. 
Even though 6 (r) is concentrated in the surface region, the average density seen through 
this reaction still amounts to 0.069 fm~^ (0.080 fm~^) or 32% (45%) of the central density 
in the case of lsi/2 knockout from ^^C (^'^Ca). In the case of emission from the Fermi level, 
on the other hand, the average density is only 12% (13%) of the central density for a ^^C 
(''°Ca) target. 

Also, the IFSI factors for neutron emission are almost identical to the corresponding IFSI 
factors for proton knockout and, as expected, the overall effect of IFSI is more pronounced 
for heavier target nuclei. 

E. Comparison between ROMEA and RMSGA calculations 

In this subsection, we investigate the sensitivity of the computed IFSI factors to the 
adopted parametrizations for the optical potentials and compare the ROMEA results with 
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FIG. 7: Contribution to the 2p) cross section b (r) as a function of r. The upper figures present 
the results obtained after setting the IFSI factor Sipsiif) equal to one in Eq. HB4|) . whereas in 
the lower panels the ROMEA calculations using the EDADl optical potential are depicted. The 
dashed (dot-dashed) curves show the result for emission from the Fermi (lowest-lying lsi/2) level. 
The baryon density p(r) is also shown (solid curve). The ordinate is given for p{r\ The 5 (r) 
are plotted in units of fm^ up to an arbitrary scaling factor. The kinematics was Tpi = 1 GeV, 
Tfci = 870 MeV, Qx = 13.4°, and ^2 = 67°. 



23 



the RMSGA predictions. As can be seen in Fig. |Hl the IFSI factor depends on whether A- 
dependent (EDAD1/EDAD2) or A-independent (EDAI) fits for the potentials are selected, 
but the global features are comparable. Figure El as contrasted to Fig. |21 demonstrates that 
the RMSGA method adequately describes the ISI of the incoming proton and the FSI of 
the scattered proton. However, the discrepancies between ROMEA and RMSGA become 
significant in the calculation of the FSI of the ejected nucleon (note the different scales in 
the bottom left panels of Figs. Olandini), and therefore also in the complete IFSI factor. 
The noted difference is attributed to the low ejectile kinetic energy {Tk2 ~ 114 MeV for 
the specific case of Fig. IHl and comparable values for knockout from other levels and other 
nuclei). At such low energies, the RMSGA predictions are not realistic due to the underlying 
approximations, mostly the postulation of linear trajectories and frozen spectator nucleons. 
So, for the kinematics discussed here, the ROMEA method is to be preferred over the 
RMSGA one, as the latter overestimates the distortion for the low-energetic ejectile nucleon. 



IV. NUMERICAL RESULTS FOR A{p,pN) DIFFERENTIAL CROSS SECTIONS 
A. The PNPI experiment 

The PNPI experiment 



42] was carried out with an incident proton beam of energy 1 GeV. 
The scattered proton was detected at 6i = 13.4° with a kinetic energy between 800 and 
950 MeV, while the knocked-out nucleon was observed at 62 = 67° having a kinetic energy 
below 200 MeV. 

Figures ITUHT^ display a selection of differential cross section results as a function of the 
cinetic energy of the most energetic nucleon in the final state. The EDAI optical potential 
3ll | was used for the ROMEA calculations. The other parametrizations of Ref. produce 
similar predictions, whereas the RMSGA approach fails to give an adequate description of 
the data because of the low kinetic energy of the ejected nucleon. Since the experiment of 
Ref. j3l only measured relative cross sections, the ROMEA results were normalized to the 
experimental data. 

The ROMEA calculations reproduce the shapes of the measured differential cross sections. 
Furthermore, comparison of the RPWIA and ROMEA calculations shows that the effect of 
the IFSI is twofold. First, IFSI result in a reduction of the RPWIA cross section that is 
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FIG. 8: The sensitivity of the real part of the complete IFSI factor for neutron knockout from 
^^C to the adopted choice for the parametrization of the optical potentials. Results of ROMEA 
calculations with the EDADl (solid curve), EDAD2 (dashed curve), and EDAI (dot-dashed curve) 
optical potentials are shown. Kinematics as in Fig. El 
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FIG. 9: (Color online) As in Fig. El but now using the RMSGA method. 



both level and A dependent. From the figures it is clear that ejection of a nucleon from 
a deeper lying level leads to stronger initial- and final-state distortions. This refiects the 
fact that the incoming and outgoing nucleons encounter more obstacles when a deeper lying 
bound nucleon is probed. The A dependence also conforms with our expectations, i.e., the 
IFSI effects are larger for heavier nuclei. Besides the attenuation, the IFSI also make the 
measured missing momentum different from the initial momentum of the struck nucleon. 
As can be inferred from Fig. ^1 this momentum shift leads to an asymmetry between the 
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FIG. 10: Differential cross section for the ^^C(p, 2p) reaction. The sohd curve represents the 
ROMEA calculation, whereas the dashed curve is the plane wave result reduced by the indicated 
factor. The ROMEA results are normalized to the data. Data points are from Ref. [4^. The 
magnitude of the experimental error bars is estimated to be of the order of 5-10%. 




FIG. 11: As in Fig. ^1 but for the ^^0{p,pn) reaction. 



positive and negative missing-momentum side of the momentum distribution. Note that 
positive missing momentum corresponds to = ki sin6'i + /c2 sin ^2 cos 02 > 0. 

B. The TRIUMF ^He(p, 2p) experiment 

Finally, we present some results for the ^He(p, 2p) reaction at an incident proton energy 
of 250 MeV. Figure El compares the data from the TRIUMF experiment with ROMEA 
calculations using the optical potential of van Oers et al. P| . The typical shape for knock- 
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FIG. 12: As in Fig. ^Jbut for the ^'^Ca(p,pn) reaction. 
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FIG. 13: The ^^C(p, 2p) momentum distribution for the 1^3/2 state as a function of the missing 
momentum. The sohd (dashed) curve are ROMEA (RPWIA) calculations. 



out of an s-state proton is reproduced by the ROMEA predictions. This fair agreement 
of the ROMEA results with the data demonstrates that our ROMEA model also works 
satisfactorily at lower incident energies. 



V. CONCLUSIONS 



A relativistic and cross-section factorized framework to describe the IFSI in quasielastic 
A{p,pN) reactions has been outlined. The model, which relies on the eikonal approximation, 
can use either optical potentials or Glauber multiple-scattering theory to deal with IFSI. 
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FIG. 14: Differential cross section for the ^^116(^,2^) reaction at the angle pairs (40°, 40°) and 
(45°, 35°) at 250 MeV. The sohd (dashed) curves refer to ROMEA (RPWIA) results. The data are 
from Ref. |9||. 



Thanks to the freedom of choice between these two substantially different techniques, our 
model is expected to be applicable at both intermediate and high incident energies. 

The properties of the IFSI factor have been investigated for an incident proton energy of 
1 GeV. Not surprisingly, the strongest attenuation occurs in the nuclear interior and heavier 
target nuclei are found to induce larger IFSI effects. Also, the surface-peaked character of 
the A{p,pN) reaction was clearly established to be a consequence of the IFSI. Whereas the 



different types of optical-potential sets contained in Ref. [31| yield comparable IFSI factors, 
the RMSGA calculations exhibit an unrealistic behavior for these kinematics. 

The ROMEA model has been used to calculate cross sections for the kinematics of two 
different experiments : quasielastic proton scattering from ^^C, ^^O, and ^°Ca at 1 GeV, 
and ^He(p, 2p) scattering at 250 MeV. The predictions are shown to reproduce the shape of 
the data reasonably well at both incoming energies, thereby providing support for the wide 
applicability range of our model. 
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Although the RMSGA approach was deemed unsuitable for the kinematics discussed 
here, it should prove useful when trying to describe nuclear transparency data. Work in this 
direction is in progress. 
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APPENDIX A: RELATIVISTIC BOUND-STATE WAVE FUNCTIONS 



For spherically symmetric potentials, the solutions 0^ (0 to a single-particle Dirac equa- 



tion have the form 



r, cr 



r, cr 



r 



(Al) 



where n denotes the principal, k and m the generalized angular momentum quantum num- 
bers. The y±Km are the spin spherical harmonics and determine the angular and spin parts 
of the wave function. 



minis 



y-Kmi^,^) = Yl yiT^i^^s I jm^ lj^^(fi)xi„^(a) 



minis 



(A2) 



with 



J = IN 



I 



K>0 



1 = 23-1 



K — 1, K > 



(A3) 



2 I -(k-M), K<0 I -K, K<0 

The Fourier transform of the relativistic bound-nucleon wave function is given by 

\l ,'9^^3/2 dnnip) 3^Km(^p) 

~Si^ fnnip) y-KmiS^p) 

with S'k = k/ The radial functions and in momentum space are obtained from 
their counterparts in coordinate space : 



<t>a{p)= / dre-'^^Uc (r1 = {-1)' (27r)= 



(A4) 



gnniP) = M/ - / r dr ji{pr) , 

V vr Jo r 



(A5a) 
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fnnip) = r'dr ^^jj{pr) , (A5b) 

V TT ./n r 



with ji{pr) the spherical Bessel function of the first kind. 



APPENDIX B: RADIAL CONTRIBUTION TO THE A{p,pN) CROSS SECTION 



The differential A{p,pN) cross section ()35p is proportional to the distorted momentum 
distribution of Eq. (j^Sl)- When approximating the completeness relation (P3j) as 



i.e., neglecting the negative-energy term as in Sec. IllBl this amounts to 



dEkidVLidVL'2 



oc 



Thus, with D (r) defined as 



D{r) = j dn e"*P™-^>„, (f) Sjpsi (r1 



the function 



-I POO POO 



dr D{r)+ drD{r)){ dr D (r) + / c/r D (r; 

Jri+AR ^ ^Jo Jri+AR 



fBll 



(B2) 



(B3) 



(B4) 



represents the contribution of an infinitesimal interval in r around ri to the A{p,pN) cross 
section. This procedure also enables us to estimate the average density seen through this 
reaction as 



P 



(B5) 
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